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Abstract. High precision photometry and spectroscopy of low-mass stars reveal a variety 
of properties standard stellar evolution cannot predict. Rotation, an essential ingredient of 
stellar evolution, is a step towards resolving the discrepancy between model predictions and 
observations. 

The first rotating stellar model, continuously tracing a low-mass star from the pre-main se- 
quence onto the horizontal branch, is presented. The predicted luminosity functions of stars 
of globular clusters and surface rotation velocities on the horizontal branch are discussed. 
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1. Introduction 

Rotation has to be included in realistic stel- 
lar models. Since recent high precision pho- 
tometry and spectroscopy have become avail- 
able, the effects of slower rotation on low-mass 
stars has become detectable in the colour- 
magnitude diagram (CMD) and Luminosity 
Functions (LFs). To predict them precisely, 
low-mass stellar models need to take rotation 
into account. 



1.1. Stellar Structure 

The version of the stella r evolution code 
FRAN EC, described in iPietrinferni et al. I 

(2004), is used for our computations. 
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Rotation was included in the stellar struc- 
ture equations following the prescription of 
Kippenhahn & Thomas! dl970h . The required 
parameters /p and fj were derived analytically 
from the effective gravitational potential of the 
star approximating its equipotential surfaces 
as rotational ellipsoids. 

We have eliminated the subatmospheric re- 
gion used in the original code and increased the 
density of the outermost meshpoints in the stel- 
lar model. 



1.2. Transport Mechanisms 

Chemical elements and angular momentum 
are transported within the star during its evo- 
lution. In a convective zone (CZ) the trans- 
port is assumed to be instantaneous. The 
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abundance of chemical elements is averaged 
in a CZ and the angular momentum re- 
distributed to achieve either specific angu- 
lar momentum conservation or solid body ro- 
tation. Both di stributions are limiting case s 
introduced by ISweigart & Mengel I 1 19791) . 
Solid body rotation is chosen for CZ un- 
til the the star reaches the turn off (TO), 
after which the specific angular momentum 
is conserved (se e Palaci os & Brun 1 120061 and 
ISills & Pinsonneauitll2000h . 

In a radiative region the transport is a diffu- 
sive process, since Palacios et al. (2006) have 
shown that advective transport via meridional 
currents is less efficient than via diffusion re- 
sulting from the shear instability. This instabil - 
ity dominates rotational mixing dZahn II 19921) . 
The diffusion coeffici e nt D y is taken from 
iDenissenkov & Weiss I (120041) introducing an 
efficiency parameter F. Chemical elements are 
transported by a combination of atomic diffu- 
sion and rotational mixing. Angular momen- 
tum is diffused by rotational mixing alone. 

1.3. Braking 

The total amount of angular momentum of the 
star is changed by several processes. On the 
pre-main sequence (PMS) magnetic fields con- 
nect the star with its circumstellar disk. The 
disk allows angular momentum to be diffused 
beyond the influence of the magnetic field of 
the star, preventing t he PMS spin-up. This pro- 
cess of disk braking dShu et al. II 1 9941) is active 
until a time tdb- 

On the main sequence (MS) the stel- 
lar magnetic field increases the amount 
of angular momentum lost by its stel- 
lar wind. The angular momentum lost by 

this magnetic-wind braking is giv en by 

IChaboyer. Demarque. & Pinsonneault I (Il995l) . 
We include a function determining the effi- 
ciency to generate a magnetic field B e g and 
parametrise it by the mass of the convective en- 
velope. The shape of B e ff is determined by two 
parameters (M W ;B W ). The overall efficiency of 
the braking process is determined by /&. 

Angular momentum lost via mass loss is al- 
ready included in the code, and is dominant on 
the red giant branch (RGB). 



2. Calibration 

The modified stellar evolution code contains a 
number of free parameters. First, the two pa- 
rameters describing the initial conditions of the 
star have to be determined: the initial rotation 
rate coq and the disk braking time tb- Four ad- 
ditional parameters to be calibrated are: the 
mixing length a, the braking efficiency /k, the 
function describing the efficiency to generate 
the magnetic field B e ff, and the rotational mix- 
ing efficiency F. 

2.1. Initial Conditions 

To determine ujq we used the rotation period of 
a sample of ~1 Myr old stars from thr ee open 
cluster: IC 348 (Littlefair et a lTI l2005h . Orion 
Nebula Cluster (Herbst et al. 200lTand_the 
Orion Nebula Cluster Flanking fields (iRebull I 
2001). We limited the mass range to 0.5 M Q < 
M < 1 .3 M to later reproduce the surface rota- 
tion velocity for the Hyades. The mean value of 
the rotation period at 1 Myr is Pi Myr = (5 ±4) d. 
Appropriate values of uq for our stellar mod- 
els were chosen to reproduce P\ M yr - This wq 
is derived from stars with solar metallicity. To 
determine wq at lower metallicites, the total an- 
gular momentum of stars of the same mass, but 
different metallicity, is assumed to be the same 
at the start of their life on the PMS. 

The disk braking time is set to tdb = 
5 Myr, at which 50 % of the stars in clus- 
ters have no longer a circums tellar disk 
dBouvier. Forestini. & Allain Ill997b . 



2.2. Free Parameters 

To calibrate the mixing length (a) and /j we 
created a theoretical solar model. The pre- 
dicted rate of angular momentum loss ( ^ J for 
the model was used to confirm the adopted 
braking mechanisms. In Fig.Q]we give the rate 
of angular momentum loss for three different 
solar models with initial rotation rates con- 
sistent with observed -Pi Myr (PiMyr=5d grey 
line, PiMyr=12d thick light grey line, and 
P\ Myr= 1 d thin black line). Each of these mod- 
els reaches \M\ observed for the Sun and 
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Fig. 1. Rate of angular momentum loss (^) by 
magnetic wind braking plotted against surface rota- 
tion rate of three different solar models (P[ M yr -5 d 
grey line, PiMyr=12d thick light grey line, and 



1 Myr' 



: 1 d thin black line). 



Observed ^f values for 

dr 



the Sun (log — =0) and a theoretically derived up- 
per limit (dotted line) for the MS is shown (see text 
for details). 
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Fig. 2. The observ ed (vsini) taken from 
ISoderblom et al. I dl993l) as a function of T cr r. 
The profile of the B e g functions are given in the 
inset to the upper right, including the point (M W ;B W ) 
that defines its shape. Isochrones for (vsini) are 
shown using different B eff functions: (0.005;0.500) 
solid, (0.007;0.500) long dashed, and (0.009;0.508) 
short-dashed line. 



shown in Fig. Q] Additionally, they lie below 
a theoretical up per limit determined by a 2D 
IWeber & Davis ( 1967) equatorial wind model. 

Z? e ff was determined by predicting sur- 
face rotation velocities for the Hyades 
and comparing th em w ith observations by 
ISoderblom et al. I dl993l) . Figure |2] shows 
three predicted rotation distributions for 
different choices of the B e ff function ( 
(0.005;0.500) solid, (0.007;0.500) long 
dashed, and (0.009;0.508) short-dashed line). 
The profile of the B e ff functions are given in 
the inset to the upper right, including the point 
(M W ;B W ) that defines its shape. The best fit 
between data and models was reached with 
(0.005;0.500). 

Spectroscopic observations that in the 
past were explained in terms of efficient ro- 
tational mixing are now being challenged 
by new data and theoretical developments 
(see Spite plateau data of Asplund et al. 



Ventura and E. Carretta in this proceedings, 
for the issue of extra mixing on the RGB). 
Therefore, the rotational mixing efficiency F is 
set to an arbitrary low value of 10~ 6 , allowing 
for a marginal transport without creating any 
significant changes to the chemical abundances 
in the stellar model. 



3. Application to GC 

Using our rotating stellar models we present 
a comprehensive analysis of rotational ef- 
fects on the evolution of low-mass stars in 
lobular clusters. Past work was explora tory 



2006 presented in Bonifacio et al. 2007; 
Egglet on. Dearborn. & Lattanziol (120071) . P. 



Vande nberg. Larson. & D e Propris | | l998l) o r 
targeted specific issues (Palacios et al. 2006). 
We apply our models to typical globular cluster 
stars (Z=10~ 3 and Y=0.246) and correspond- 
ing isochrones of 10.5 and 12.0 Gyr. 



3.1. Surface Velocities 

Our isochrones predict a surface velocities 
on the SGB v sur f = 1 - 5.5kms _1 , consis- 
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tent with the observed \w ~ 3kms~' by 
iLucatello & Gratton I d2003i More precise ob- 
servations of the rotation along the SGB are 
necessary to further constrain the initial con- 
ditions and the angular momentum transport in 
CZs. 

Rotational velocities on the HB were 
determined for a 0.8 M star experiencing 
no mass loss with an optimised version of 
our rotational code that is able to evolve 
through the helium flash. We predicted v sul f = 
6.85 kms~' on the zero-age HB (ZAHB) that 
is within the range observe d for red HB stars 
dSills & Pinsonneaultll2000l) . 



3.2. Isoch rones 

Two isochrones for ages of 10.5 and 12.0 Gyr 
are shown in Fig. [3] with (dashed line) and 
without rotation (solid line). Rotation does not 
shift the MS or the TO, however, the brightness 
of the SGB at (B-V)=0.55 is slightly increased 
by 0.05 and 0.03 mag, for 10.5 and 12.0Gyr re- 
spectively. The RGB at My - 2 including ro- 
tation is bluer by 0.01 and 0.007 mag for 10.5 
and 12.0 Gyr respectively. 

Rotation increases the tip of the RGB 
(TRGB) luminosity by A log L/L Q = 0.08 and 
the helium core mass at the flash by AM^e - 
0.048 Mq. This large helium core leads to a 
brighter ZAHB luminosity of AM V = 0.326 
when including rotation, and will be discussed 
later. 



3.3. Luminosity Functions 

From the isochrones we then determine the LF 
for the case of r otation and no rotation using a 
ISalpeter I (1 19551) initial mass function. Figured 
shows the two cases are normalised at a My-4 
which is close to the TO and independent of the 
choice of initial mass function. The shape of 
the RGB bump is not changed by rotation, only 
its location is shifted to brighter luminosities 
by AM V = 0.1 for the 10.5 Gyr LF Simil ar to 
IVandenberg. Larson. & De Proprisl(ll998l) . we 
find that rotation reduces the SGB slope and 
increases the number count at the base of the 
RGB (by 1 1 % for the 12 Gyr LF). 
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Fig. 3. The 10.5 Gyr and 12 Gyr isochrones with- 
out rotation (as labelled) are plotted in a BV CMD 
as solid lines. The isochrones including rotation are 
shown as dashed lines. 



4. Discussion 

Given an isochrone without rotation that is 
consistent with observational constraints. The 
inclusion of rotation increases the helium core 
mass, this leads to a brighter ZAHB and dis- 
tances from HB fitting that are no longer con- 
sistent with parallax b ased empirical MS fitting 
dRecio-Blanco et al. Il2004h . 

A solution would be to reduce the rotation 
rate in the helium core on the RGB. This re- 
duction may result from either more realistic 
angular momentum trans port in CZs (see ex- 
ploratory work by Palaci os & Brun I [20061) or 
reducing ojq. 

We chose to reduce the initial rotation rate 
by adopting a different method to transpose so- 
lar metallicity rotation rates to lower metallici- 
ties. Instead of keeping the total amount of ini- 
tial angular momentum constant, a>o is defined 
so that these stars reach a rotational period of 
Pi Myr=5 d. This reduces wq by 60 % and leads 
to a SGB rotation velocity of 3 km s , consis- 
tent with observations. The ZAHB including 
rotation is only brighter by AMy =0.122 (com- 
parable to typical photometric errors for HB 
stars) and leads to ages younger by 1.5 Gyr. 
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Fig. 4. Three LFs are shown as a function of M v 
with an age of 12 Gyr without (solid line) and with 
rotation (dashed lines), as well as an age of 10.5 Gyr 
including the reduced initial rotation (dotted line). 



Therefore, an isochrone or LF with and with- 
out rotation can only be consistent with obser- 
vational constraints, if the rotational isochrone 
or LF is 1 .5 Gyr younger than without rotation. 
A consistent comparison between a 
10.5 Gyr old LF with the reduced o>q (solid 
line) and a 12.0 Gyr LF without (dashed line), 
is shown in Fig. |4] The prior determined 1 1 % 
increase in number count is reduced to 5 %, 
but it is still detectable. 
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